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ABSTRACT 



We present medium-resolution optical spectra of the secondary star in the high Galactic 
latitude black hole X-ray binary XTE J1118+480 and determine the abundance of Mg, Al, Ca, 
Fe, and Ni in its atmosphere. For all the elements investigated we find supersolar abundances; 
thus, we reject the hypothesis that the black hole came from the direct collapse of an ancient 
massive halo star. The compact primary most likely formed in a supernova event of a massive star 
whose nucleosynthetic products polluted the secondary star. The observed element abundances 
and their ratios can be explained using a variety of supernova models with a wide range of 
metallicities. While an explosive origin in the Galactic halo or thick disk cannot be discarded, a 
metal-rich progenitor is clearly favored by the observed abundance pattern. This suggests that 
the black hole was produced in the Galactic thin disk with a violent natal kick, propelling the 
X-ray binary to its current location and orbit. 

Subject headings: black hole physics — stars: abundances — stars: evolution — stars: individual 
(XTE J1118+480) — supernovae: general — X-rays: binaries 



1. Introduction 

The low-mass X-ray binary XTE J1118+480 
was discovered with the all-sky monitor aboard the 
Rossi X-ray Timing Explorer on UT 2000 March 
29 (Remillard et al. 2000). Throughout its out- 
burst the source remained in the low/hard state, 
one of the characteristic spectral states of an ac- 
creting black hole binary (McClintock et al. 2003). 
The system consists of a black hole with a mass es- 
timated in the range 6—8 M Q and a late-type sec- 
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ondary star of 0.1-0.5 M (Wagner et al. 2001). 

The extraordinarily high Galactic latitude (b « 
62.3°), together with its distance of 1.85±0.36 kpc 
(Wagner et al. 2001), places the system at a height 
of ~1.6 kpc above the Galactic plane. In addition, 
an accurate measurement of its proper motion cou- 
pled with its distance provides space- velocity com- 
ponents U, V which seem consistent with those 
of some old halo globular clusters (Mirabel et al. 
2001). If the system formed in the Galactic halo, 
the black hole could be either the remnant of a su- 
pernova in the very early Galaxy or the result of 
a direct collapse of an ancient massive star. How- 
ever, the galactocentric orbit crossed the Galactic 
plane many times in the past, and an alternative 
possibility is that the system formed in the Galac- 
tic disk and was launched into its present orbit as 
a consequence of the "kick" acquired in the super- 
nova explosion of a massive star (Gualandris et al. 
2005). The metallicity of the secondary star may 
provide a key to distinguish among these two pos- 
sible birth places, giving important clues to the 
formation of the black hole and the properties of 



the supernova explosion, such as symmetry, re- 
leased energy, and characteristics of the ejected 
matter (Israelian et al. 1999; Podsiadlowski et al. 
2002). 

2. Observations 

We obtained 74 medium-resolution spectra 
(X/SX w 6,000) of the secondary star in XTE 
J1118+480, in quiescence, on UT 14 February 
2004, using the 10-m Keck II telescope, equipped 
with the Echellette Spectrograph and Imager 
(ESI; Sheinis et al. 2002). The exposure time 
was fixed at 300 s to minimize the effects of or- 
bital smearing which, for the orbital parameters of 
XTE J1118+480, is in the range 0.6-26.6 kms" 1 , 
smaller than the instrumental resolution of ~ 50 
kms -1 . Each individual spectrum was corrected 
for the radial velocity of the star, and the spectra 
were combined in order to improve the signal-to- 
noise ratio. After binning in wavelength in steps 
of 0.3 A, the final spectrum had an average signal- 
to-noise ratio of 80 in the continuum. The data 
cover the spectral range 4000-9000 A and clearly 
show the characteristic emission lines of accreting 
low-mass X-ray binaries (Balmer series and Ca II 
near-infrared triplet, He I 5876 A, He I 6678 A), 
superimposed on the typical photosphcric spec- 
trum of a late-type star. 

3. Chemical Analysis 

A comparison of the observed spectrum with 
the spectra of ten template stars (K0V-M2V) ob- 
tained with the same instrument allows us to clas- 
sify the secondary as a mid to late K-type star. 
In Figure 1 we display two spectral regions con- 
taining some relevant lines for our analysis as well 
spectrum of a template star of similar spec- 
tral type for comparison with available detailed 
chemical analysis ( Allende Prieto et al. 2004) . We 
also show synthetic line profiles for different stel- 
lar abundances computed with the local thermo- 
dynamic equilibrium (LTE) code MOOG (Sneden 
1973), adopting the atomic line data from the 
Vienna Atomic Line Database (VALD; Piskunov 
1995) and using a grid of LTE model atmospheres 
(Kurucz 1993). 

In order to perform the chemical analysis of 
the secondary star we used a technique which 
combines a grid of synthetic spectra and a %2- 



minimization procedure that includes Monte- 
Carlo simulations (Gonzalez Hernandez et al. 
2004, 2005). First, we inspected the observed 
spectrum in order to select the most suitable fea- 
tures for a chemical abundance determination. 
We identified nine spectral features containing in 
total 30 lines of Fe I and 8 lines of Ca I with ex- 
citation potentials between 1 eV and 5 eV. The 
oscillator strengths of the relevant spectral lines 
were checked via spectral synthesis against the 
solar atlas (Kurucz et al. 1984). We then gener- 
ated a grid of about one million synthetic spectra 
of each of these features, varying as free parame- 
ters the star effective temperature (T c ff), surface 
gravity (logg), and metallicity ([Fe/H]), together 
with the veiling from the accretion disk which was 
assumed to be a linear function of wavelength, 
and thus described by two additional parame- 
ters. Iron abundances were varied in the range 
— 1.5 < [Fe/H] < 1 whereas the Ca abundance 
was fixed, for each given iron abundance, accord- 
ing to the Galactic trend of Ca (Bensby et al. 2005) 
for [Fe/H] < 0, and fixed to [Ca/Fe] = in the 
range < [Fe/H] < 1. A rotational broadening 
of 100 kms -1 and a limb darkening e = 0.8 were 
adopted. The microturbulcnce (£) was computed 
using an experimental expression as a function of 
effective temperature and surface gravity (Allende 
Prieto et al. 2004). 

We compared, using a x2-minimization pro- 
cedure, this grid with 1000 realizations of the 
observed spectrum. Using a bootstrap Monte- 
Carlo method, we found the most likely values 
T cff = 4700 ± 100 K, log( 5 /cm s 2 ) = 4.6 ± 0.3, 
[Fe/H] = 0.2 ± 0.2, and a disk veiling (defined 
as -Fdisk/-Ptotai) of less than 40% at 5000 A and 
decreasing toward longer wavelengths. The (ler) 
uncertainty in the iron-abundance determination 
takes into account the uncertainties in the stellar 
and veiling parameters. The effective temperature 
and surface gravity are consistent with previous 
spectral classifications and similarly the reported 
veiling values (Torres et al. 2004). Using the de- 
rived stellar and veiling parameters, we analyzed 
several spectral regions where we had identified 
various lines of Fe, Ca, Al, Mg and Ni. Abun- 
dances of all the elements are listed in Table 1. 
The la uncertainty in the abundance determina- 
tion takes into account uncertainties in the stellar 
and veiling parameters. 
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Remarkably, we find a metallicity higher than 
solar, which is extremely atypical of halo stars (Al- 
lende Prieto et al. 2006). In Figure 1 we show 
the best-fit synthetic spectrum to various features 
in two different spectral regions, and for compar- 
ison, a model with twenty-five times lower metal 
content. Notice the inadequacy of low-mctallicity 
models to reproduce the observed features, even 
in the extreme case that no veiling is considered 
(adding veiling would make the discrepancy much 
worse). An iron abundance of [Fc/H] = —1.2 
is more than 6a away from the best-fit solution 
([Fe/H] = 0.18), and hence very unlikely. We 
have also found that abundances of Al, Ca, Mg, 
and Ni arc higher than solar (see Tabic 1). In 
Figure 2, we show that the abundance ratios of 
these elements with respect to iron are consistent 
with those of stars in the solar neighborhood from 
Gilli et al. (2006). We have also determined an 
upper limit to the Li abundance loge(Li)LTE = 
log[AT(Li)/Ar(H)] L TE + 12 < 1.61 ± 0.25 using the 
Li 6708 A line. This value seems to be lower 
that typical high Li abundances measured in other 
late-type secondary stars in soft X-ray transients 
whose origin is still an open question (Martin et 
al. 1994). 

4. Discussion and Conclusions 

If we include the metallicity distribution of halo 
and thick-disk stars (Allcnde Prieto et al. 2006) 
and thin-disk stars (Allende Prieto et al. 2004) in 
the equations based on kinematics for establish- 
ing the relative likehoods of belonging to the halo, 
thick disk, or thin disk (Bensby et al. 2003), the 
probability that a star with the Galactic space ve- 
locity components of this system (Mirabel et al. 
2001; U = -105 ± 16 kms" 1 , V = —98 ± 16 
kms -1 , W = —21 ± 10 kins" 1 ) and metallicity 
[Fe/H] = 0.18 belongs to the Galactic halo is less 
than 0.1%. Moreover, the kinematics alone sug- 
gest thick-disk rather than halo membership, al- 
though the high metallicity of the secondary star 
favors thin-disk membership. If, however, the pro- 
genitor was a massive star in the halo or thick 
disk, the high metallicity of the secondary rules 
out the hypothesis that the black hole was formed 
by direct collapse of the massive star; instead, a 
supernova explosion origin is strongly suggested. 

Given its present orbital distance from the black 



hole (a c 3 R©, Wagner et al. 2001), it is plausi- 
ble that the secondary star captured a significant 
fraction of the matter ejected in a supernova explo- 
sion. The chemical composition of the secondary 
may provide crucial information on nucleosynthe- 
sis in the progenitor and the formation mechanism 
of the black hole. We consider two possible scenar- 
ios for the origin of the compact object: either it 
formed in the Galactic halo or thick disk as a result 
of the explosion of a metal-poor massive progen- 
itor which enriched the secondary star from the 
typical abundances of halo/thick-disk stars up to 
the observed supersolar values, or alternatively, it 
formed in the Galactic thin disk with a natal kick 
imparted during the supernova explosion which 
propelled the binary into its current orbit. 

In the first case, the similarity with the kine- 
matics of halo and thick-disk stars makes unnec- 
essary a significant kick during the black hole for- 
mation process. However, since the typical metal- 
licities of halo stars and thick-disk stars are sig- 
nificantly lower than solar, it is required that the 
secondary captured enough matter from the ejecta 
to reach the current abundances. A metal-poor ~1 
M© secondary star, initially placed at an orbital 
distance of ^6 R© (after tidal circularization of 
the orbit), would need to capture roughly 5-10% of 
the matter ejected in a spherically symmetric core- 
collapse supernova explosion of a 16 M© helium 
core to achieve the observed iron abundance. We 
have considered supernova and hypernova models 
of metal-poor progenitors with different masses, 
metallicities, mass cuts (i.e., the mass above which 
the matter is expelled at the time of the supernova 
explosion), fallback (i.e., amount of mass which 
is eventually accreted by the compact core), and 
mixing (Umeda & Nomoto, 2002, 2005; Tominaga, 
Umeda & Nomoto 2006, in preparation), and as- 
sumed that all the fallback material is well mixed 
with the ejecta. We find that a sufficiently large Fe 
enrichment is possible for mass cuts in the range 
2-4 M© [and explosion energies (1 - 30) x 10 51 
ergs] . The abundance ratios of the other elements 
are only marginally reproduced by these models. 
An origin of XTE Jl 118+480 in the Galactic halo 
or in the Galactic thick disk cannot be discarded 
by the present observations, nor can it be con- 
firmed. Further higher-quality observations and 
extensive model calculations are required to fully 
explore this possibility. Further details of analy- 
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sis will put together in Gonalez Hernandez et al. 
(2006, in preparation). 

In the second scenario, the system had to ac- 
quire a peculiar space velocity of ~ 180 km s" 1 , to 
change from a Galactic thin disk orbit to the cur- 
rently observed orbit (Gualandris et al. 2005), re- 
quiring an asymmetric kick. Such kicks imparted 
during the birth of nascent neutron stars, due to 
asymmetric mass ejection and/or an asymmetry in 
the neutrino emission (Lai et al. 2001), have been 
proposed to explain the large transverse motions 
of neutron stars in the plane of the sky (Lyne & 
Lorimer 1994). The black hole could have formed 
in a two-stage process where the initial collapse 
led to the formation of a neutron star accompa- 
nied by a substantial kick and the final mass of the 
compact remnant was achieved by matter that fell 
back after the initial collapse as proposed to ex- 
plain the origin of the black hole in Nova Sco 1994 
(Podsiadloski et al. 2002; Brandt et al. 1995). 

Spherically symmetric explosion models are 
able to explain the observed metal enrichment. 
A modest amount of ejecta (~ 0.02) M Q captured 
by a secondary with initial solar abundance is suf- 
ficient if vigorous mixing (Kifonidis et al. 2000) 
between the fallback matter (which is necessarily 
large given the mass of the black hole) and the 
ejecta took place. In this scenario the black hole 
would be formed in a mild explosion with fallback 
as in collapsar models (MacFadyen et al. 2001) 
associated with rapidly rotating massive stars. 
However, in the pure spherical supernova explo- 
sion scenario for XTE J1118+480, the maximum 
allowed ejected mass to keep the secondary star 
gravitationally bound is AM « 8 M© . This leads 
to a system velocity of ^60 kms -1 (there could 
be extreme cases with system velocities as high 
as ~ 100 kms -1 requiring very special model pa- 
rameters such as higher secondary and black hole 
masses just after the explosion; see further details 
in Sec. 3.2 of Gualandris et al. 2005), much lower 
than the observed value, and a neutrino-induced 
kick is required if we are to explain the kinematics 
of the system. 

The chemical composition of the ejecta in a 
non-spherically symmetric supernova explosion is 
strongly dependent on direction. In particular, if 
we assume that the jet is collimated perpendicular 
to the orbital plane of the binary, where the sec- 
ondary star is located, elements such as Ti, Ni, and 



Fe are mainly ejected in the jet direction, while Al, 
O, Si, S, and Mg are preferentially ejected near the 
equatorial plane of the helium star (Maeda et al. 
2002). Using predictions for an asphcrical explo- 
sion model of a 16 M He core metal-rich progen- 
itor, we can explain the observed abundances in 
the secondary. Complete lateral mixing (Podsiad- 
lowski ct al. 2002) is required to account for the 
similar enhancement of Mg and Ni and the ob- 
served abundances can be reproduced for all mass 
cuts in the range 2-8 M . 

It is therefore plausible that the black hole 
in XTE J1118+480 formed in the Galactic thin 
disk from a massive metal-rich progenitor and was 
launched into its current orbit either by mass ejec- 
tion in an asymmetric supernova/hypernova ex- 
plosion or by a neutrino-induced kick. In addi- 
tion, ultraviolet observations of the accretion disk 
in this system suggest that the material accreted 
onto the compact object is substantially CNO pro- 
cessed (Haswell et al. 2002), indicating that the 
zero-age mass of the secondary star could have 
been ~ 1.5 M Q . By means of binary evolution cal- 
culations this may constrain the age of the system 
in the range 2-5.5 Gyr (Gualandris et al. 2005). 
Future chemical studies of the secondary star dur- 
ing periods of quiescence may provide accurate 
abundances of elements (e.g., C, Ti, Si) whose lines 
either are not available in the observed spectrum 
or are present in spectral regions where the signal- 
to-noise ratio is too low for an accurate chemical 
analysis. These may reveal further details of the 
formation mechanism of the black hole in this sys- 
tem. 
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Fig. 1. — Best synthetic spectral fits to the 
ESI spectrum of the secondary star in the 
XTE JII18+480 system (second and bottom pan- 
els) and the same for a template star (prop- 
erly broadened) taken from Allcnde Prieto et 
al. (2004) shown for comparison (top and third 
panels). Synthetic spectra are computed for 
typical abundances for a halo star ([Fe/H] = 
— 1.2, dashed-dotted blue line), solar abundances 
([Fe/H] = 0, dashed green line), and best-fit abun- 
dances (solid red line). In addition, note that 
for solar and best-fit abundances we have applied 
to the synthetic spectra the corresponding val- 
ues for the veiling according to the solution found 
with the fitting procedure. However, for the low- 
metallicity synthetic spectra, we have not assumed 
any veiling. 
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Fig. 2. — Abundance ratios of the secondary star 
in XTE J1118+480 (blue wide cross) in compar- 
ison with the abundances of G and K metal-rich 
dwarf stars. Galactic trends were taken from Gilli 
et al. (2006). The size of the cross indicates the 
uncertainty. Filled and empty circles correspond 
to abundances for planet host stars and stars with- 
out known planet companions, respectively. The 
dashed-dotted lines indicate solar abundance val- 
ues. 
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Table 1 

Abundances in the secondary star 



Parameter 


Al 


Ca 


Mg 


Fo 


Ni 


[X/H] 
A[X/H] 


0.60 
0.20 


0.15 
0.23 


0.35 
0.25 


0.18 
0.17 


0.30 
0.21 



Note. — Element abundances ( calculated assum- 
ing LTE) arc [X/H] = log[JV(X)/W(H)] Bta r - 
log[AT(X)/Ar(H)] S un, where N{X) is the number 
density of atoms. Uncertainties, A[X/H], arc ler and 
take into account uncertainties in the stellar and veiling 
parameters. 
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